Two-scalar-field model for the interaction of dark energy and dark matter 
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In this paper, we study the effects of an interaction between dark matter and dark energy through 
a two scalar field model with a potential V(4>, x) = e_A ^-P(0>x)i where P((f>,x) is a polynomial. 
We show that the cosmic expansion dynamics of the Universe is reproduced for a large range of 
the bare mass of the dark matter field and that there exist solutions with transient accelerated 
expansion. A modification in the exponential behavior of the potential is studied, with important 
physical implications, including the possibility of more realistic transient acceleration solutions. 
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I. INTRODUCTION 

A large number of models has been proposed to ex- 
plain the dark sector of the Universe (see Refs. [l|, for 
reviews on dark energy (DE) and dark matter (DM) , re- 
spectively). Most of those models assume that the dark 
components are noninteracting and treat them as fluids. 
However, there are neither theoretical arguments forbid- 
ding such an interaction, nor there exist sufficient ob- 
servational results to rule it out. It is then just natural 
to study the general situation in which dark matter and 
dark energy are coupled, in order to gain a deeper insight 
into the nature of these componentflj. This is motivated 
by the theoretically appealing idea that the full dark sec- 
tor can be treated in a single framework. Moreover, the 
fact that the present energy densities of dark energy and 
matter are observed to be the same order of magnitude, 
suggests a connection between them. 

An intimate connection between DE and DM is natu- 
rally expected in unification models, such as for instance 
the Chaplygin gas model and its generalizations 0-0] ■ 
Actually, in the context of this model, an assumption 
about the equation of state (EOS) of DE allows to extract 
an explicit interaction between the dark components Q . 
A map between the generalized Chaplygin gas (GCG) 
model and the interaction model to be discussed in the 
following paragraph can be found in Ref. Q . 

A general way to describe the DM-DE interaction is to 
introduce an energy exchange term Q in the conservation 
equations as follows: 



Pde + 3H(p de + p de ) = Q , p drn + 3Hp dm 



(1) 



One may phenomenologically study this interaction by 
withholding any assumptions about the nature of the 
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dark sector and treat it straightforwardly as a two- 
component fluid. The coupling Q is usually taken to 
be of the form Q = S de Hp de + S dm Hp dm , where H is the 
expansion rate and Si are coupling terms. This treatment 
is encountered in many observational studies [H, [l(| ■ 

An alternative path to study the interaction assumes 
that dark energy can be described by a scalar field 
4> in interaction with a fluid, the so-called interacting 
quintessence model: this is the case for the models of 
Refs. [ill, EH ! in which the coupling is chosen to be 
Q = f{<p)<f>p, where f(4>) is a generic function. A simi- 
lar mechanism is the so-called chameleon model In 
these cases the field interacts with every component of 
the Universe, leading to observable effects in solar sys- 
tem tests of gravity. A simple modification is possible 
using a much smaller coupling for baryons than for DM 
as in Ref. [l4j , or simply substituting p with p dm , so that 
DE cou ples only to DM. That is the case for the model 
in Ref. [15| , in which the quintessence potential and the 
interaction term are derived from scaling assumptions. 

A more fundamental approach to tackle the interaction 
treats DE and DM as fields, which abandons the need for 
fluids in the treatment of these components. Usually this 
is achieved through two new scalar fields, (f> for DE and 
X for DM [16H18| . An interaction potential Vi nt ((j), x) is 
then introduced to account for the energy exchange. Ul- 
timately these models also lead to the coupling of the 
interacting quintessence scenario, as long as the inter- 
action is in the form of a DM mass term. This is to 
be expected, as will be derived in the following sections. 
These are then similar to the so-called VAMP models [l9[ 
in which a particle is introduced whose mass varies with 
the quintessence fielc0. 

The advantages of this type of approach are manifold: 
the full set of coupled equations can be found from an 
action and consequently the functional form of the EOS 
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parameter and the DE-DM coupling is fully determined. 
Furthermore, this is an elegant and straightforward way 
to link DM and DE with more fundamental physics mod- 
els from which these components might stem. 

A relevant property of these interaction models is the 
existence of an extra force between dark matter parti- 
cles, not present in noninteracting models. This force 
can influence structure formation, creating an extra bias 
between baryon and dark matter fluctuations, which may 
in principle be measured through tests of the equivalence 
principle [9j, |2l|. Furthermore, the varying dark matter 
mass can also have an influence in the anisotropy spec- 
trum of the CMB since, among other effects, it may alter 
the ratio between dark matter and radiation densities at 
last scattering. This severely constrains the simplest in- 
teraction models, in which the DM mass is a linear func- 
tion of the DE field, but its effects are still to be studied 
in detail for more complex models (22[. 

In this paper we present a scalar field interaction model 
with an interacting potential V((j),x), that incorporates 
features of some quintessence models inspired in funda- 
mental physics theories. We then characterize the phys- 
ical solutions and ascertain the role of the interaction 
term in the cosmological evolution. The main equations 
of the model are derived in Sec. [TTJ Their numerical so- 
lutions are presented in Sec IIII1 along with the relevant 
physical results. Section [IV] concludes with a discussion 
of the obtained results and a brief outlook on future de- 
velopments. 



and isotropic, leading to the following field equations: 



• dV 
+ 3H4 + — = , 

ocp 

dV 

X + 3Hx + -jt = , 
ox 



(4) 
(5) 



where H = a/ a is the expansion rate. From the stress- 
energy tensor we obtain the usual expressions for the 
pressure and energy density, 



pd = \tf + \x 2 + V(4, x) 



(6) 



We also consider a Universe filled with perfect fluids for 
matter and radiation, which we consider to be uncoupled 
and consequently to evolve as p m oc a~ 3 and p r oc a~ 4 , 
respectively. Finally, we consider the Friedmann equa- 
tion, 



H 2 = \ ( Pm + Pr + \ft + \x 2 + V^, X) 



(7) 



We introduce the density parameters fli = pi/3H 2 , in 
terms of which the usual deceleration parameter reads 

q = ~ = I- (i + n r + 3w d n d ) , (8) 

a z 2 

where wg = Pd/ Pd is the EOS parameter for the fields. 



II. INTERACTION MODEL 



B. Interaction Potential 



A. Basic equations 



We consider two interacting canonical real scalar fields 
and x, whose Lagrangian density is given bj{3 



£d = -^{d^dA + d^xduX) ~ V^, X ) • (2) 

Guided by the cosmological principle, we assume the ge- 
ometry of the Universe to be given by a flat Robertson- 
Walker metric with line element 



ds 2 = -dt 2 + a 2 (t) (dr 2 + r 2 dttl 



(3) 



with a(t) being the scale factor, normalized so that at 
present a(to) = 1 for convenience, and dfl^ 2 the line el- 
ement for the 2-dimensional sphere S 2 . For the same 
reason, we consider both scalar fields to be homogeneous 



3 We use the (— , +, +, +) metric signature and natural units with 
h = c = 8ttG = 1 . As a consequence all masses come in terms of 



We shall be interested in studying the following poten- 
tial: 

V^,x) = ^P{^x) + \m 2 X 2 , (9) 

where P(cj>, x) is a polynomial in (f> and \ and m is 
the dark matter bare mass. Such exponential couplings 
are inspired from fundamental theories like string or M- 
theory, or N = 2 supergravity in higher dimensions [23l - 
[25| . Notice that the interaction of chiral superfields in the 
context of N=l supergravity inflationary models [2H [27j 
has many common features with the present model. Fur- 
thermore, the exponential term for dark energy is the 
simplest way to vary its contributions from very high 
energy to the present and to respect the bond for nucle- 
osynthesis [28j. Hence under these conditions, one con- 
siders the general interaction term with an exponential 
multiplied by a polynomial of <j> and x- 

The polynomial P(4>, x) can be separated into the in- 
teracting and noninteracting terms, P{<j),x) = P^i^) + 
Pint{4>iX)- F° r t ne noninteracting part, we choose the 
potential first studied in Ref. (29j . 



the reduced Planck mass, M p = Mpi/\^8n. 



P (0) = A+(0-0 o ) 2 



(10) 
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As for the interacting term, an obvious choice is to re- 
quire the x ne ld to be equivalent to a fluid of nonrela- 
tivistic matter, i.e. with negligible pressure. We recall 
that, according to Rcf. [30], scalar field oscillations in a 
potential V(x) — a X n with frequency (i.e. mass) much 
greater than the expansion rate H, behave like a fluid 
with an average EOS given by 



(Px) = 



(p* 



(11) 



Notice that this is equivalent to the virial theorem for 
power-law potentials, i.e. (5X 2 ) = § (V(x))- Thus, in 
order to ensure that x i s pressureless at all times, we 
must set n = 2. A main feature of the present model is 
that the interaction with the field <\> leads to an oscillation 
with varying frequency. Note however that, unlike pre- 
heating models (see for example Ref. 31]), which exhibit 
parametric resonance, the frequency here changes slowly. 
In this case the computations discussed in Ref. [30j for a 
constant frequency hold. To finish our discussion of the 
potential we rewrite it with an explicit DM varying mass 
term in terms of the DE field, 



with 



V(<I>,X) = Vie (</>) + v dm (4, x ) 



V dm {^x) = \M 2 {cp)x 2 , 



(12) 



(13) 



where the mass function M 2 (<fi) is given in our model by 
M 2 (0) = m 2 + 2P(0)e~ A * and the polynomial for P is 
written as 



P(</>) = B + C(f> + D(j) 2 , 



(14) 



where B, C and D are order unit parameters in terms of 
the appropriate powers of the reduced Planck mass. 



C. Average Evolution Equations 

Given the high frequency of the oscillations, it is rather 
infeasible to integrate the x equation numerically. For 
that reason we consider only averages of the field. In 
particular, we shall derive the equation for the dark mat- 
ter density and work with that instead. First, we define 
the dark matter density and pressure from the EOS found 
in Eq. (JTTJ) , 



Pdm = \± 2 + \m 2 ($)x 2 , 

Pdm = \± 2 - \M 2 ^)X 2 ■ 



(15) 



By construction, their averages over an oscillation cycle 
read 

(p dm ) = (x 2 )=M 2 (4)( X 2 ) , (p dm )=Q. (16) 



Next, we multiply Eq. © by x an d insert a term 
fiKmW t0 obtain 



Taking the average yields 



1 -dM 2 ((j)) , 
3 H P*m-^^f l (x 2 )=0 



(17) 



(18) 



where we have written (p dm ) as p dm and (pdm) as p dm , 
since the density is not sensible to the oscillations, to 
a good approximation. This can be easily seen by as- 
suming the rapid oscillations of x(t) are described by a 
sinusoidal function, and hence the density depends only 
on the amplitude of the oscillations, which is not affected 
by a cyclic average. 

Substituting the average of x 2 given by Eq. (fl"6|). we 
obtain 



3Hp dm — - 1 



1 dM 2 ((/)) 



2 r A/ 2 ((/>) d(j> 



-Pdn 



(19) 



So, as previously mentioned, the equivalence relation be- 
tween coupled quintessence and the field theory approach 
is established via the relationship 



ldlnM 2 (0) 



(20) 



Furthermore, Eq. (jl9)) can be formally solved as a func- 
tion of (p, through 



Pdm(4>, a) = n a 3 M(0) , 



(21) 



where no is an integration constant. Notice this corre- 
sponds to the statement that p = nM, with M being the 
DM mass and n the number density, proportional to a~ 3 . 
With this solution, the dynamics is reduced to a single 
differential equation for cf>: this can be obtained from Eq. 
dU, with V being replaced by an effective potential V e g 
given by 



Vos(<p, a) = V de (4>) + p dm {<i>, a) . 



(22) 



These equations are valid as long as M 2 (<fi) 3> H 2 , 
otherwise the oscillation regime is not relevant and we 
must also solve Eq. (0. 



D. Modified Potential 

Having defined the potential and derived the relevant 
equations, we are now ready to draw some general con- 
clusions about the importance of each of the terms of 
the potential. As will become clear below, such results 
motivate a modification of the potential. 

First, notice that for a sufficiently large (f> we have 
M 2 (4>) m m 2 . At this regime, examining p dm and its 



4 



derivative, we see that the interaction is irrelevant. On 
the other hand, for small values of 4>, the term with m 2 
can be neglected. Assuming the polynomial P{4>) to be of 
order one, the transition value 4> c between the two phases 
can be estimated by setting m 2 e A * = 1, which results in 



A 



In m 



(23) 



Thus, for <j) < <j) Cl the interaction is relevant, becom- 
ing subdominant as the value of the scalar field grows. 
We are interested in studying the late time behavior of 
the Universe, near the stage of accelerated expansion. 
It is relevant then to estimate whether the interaction 
is important at late times. The value of the DE field 
near the present 0(0) can be estimated by assuming that 
Pdeo ~ Vde(4>(0)) and that it is close to the critical density 
p c o, which gives 



0(0) 



hi/dco 



(24) 



Requiring that 
the bare mass, 



> 0(0), yields a rather low bound for 



(25) 



present, implying the absence of dark matter as such in 
the Universe in the past. 

Confronted with these problems, we modify the model 
so to allow for a difference in the behavior of both expo- 
nentials, i.e. we choose instead, 



Vfax) =e~ X * (A+(0-0 o ) 2 ) 



(28) 



+e-^P(0) X 2 + -m 2 X 2 



with A 7^ A. This relaxes the constraint on m to the less 
strict condition 



-<P C A /2A 



(29) 



This modification evades the problem associated to the 
onset of oscillations, as these start while the interaction 
is still relevant. 

A different modification could have been made, by dis- 
carding the assumption that the parameters B, C and 
D in P{4>) are 0(1) in terms of M p . However, a solu- 
tion to the problems mentioned above would require that 
they are increased by several orders of magnitude: this 
is rather unnatural, since they are already at the Planck 
scale. 



Thus, this analysis hints that the effects of the interac- 
tion will not be detected at the present unless the DM 
bare mass is unnaturally small. If P(<fi) is O(10 s ), this 
estimate increases by roughly 2s orders of magnitude, 
which would only shift the naturalness problem to P(0). 

Another important situation is the onset of the os- 
cillatory phase: we must establish the field value for 
which M 2 (0) > H 2 . In order to obtain it, we use the 
well-known result [l| that exponential potentials lead to 
scaling solutions. Albeit this is not strictly valid in our 
model, since it is not a pure exponential, we proceed by 
assuming a scaling behavior before falls in the minimum 
of the potential, since the polynomial P(<f>, x) varies little 
during that stage. In that case, we have 



a, 



3(w + l) 
A 2 



V de 



9(w 



A 2 



^H 2 



(26) 



in which w is the effective EOS parameter for the com- 
bination of all the components. We can then rewrite the 
oscillation condition as 



9(w + 1) m 2 e x * + 2P 
A 2 K 



> 1 • 



(27) 



We recall that if the scaling occurs during nucleosyn- 
thesis, then A > 10 [2^|. With such a value for A and 
assuming that P ~ P^, the l.h.s. is always less than 
unity when the interaction is relevant, meaning that dur- 
ing that stage the field x has not begun oscillating; con- 
versely, one expects oscillations to start as the interaction 
becomes unimportant. In particular, for the threshold 
mass of Eq. (|25p. the field may not oscillate until the 



III. NUMERICAL RESULTS 

Let us start by rewriting the equations in terms of the 
number of e-folds N — In a and, for convenience, use the 
rescaled variables of Ref. [33 , 



H 

H= JT e 



2.V 



*= — e 



2N 



X = -^e 2N 



(30) 



where Hq = 72 km/s/Mpc is the present value of the 
Hubble constant. Thus, Eqs. flD, © and (jTj) now read 



H 2 — Q m Qe N 
H(X'+X) 4 



1 -9 

+ O-o + - 
6 

e 4jv OV _ 
H 2 d<f> ~ ' 

e 4JV dV 

H 2 d X 



1 ~, 
t:X< 



a 4N 



3H 2 



V(frx) , 



(31) 



= 



where the primes denote derivatives with respect to N. 
These changes improve the numerical robustness of the 
system by shortening the range of values taken by the 
new variables. From the onset of the oscillatory phase we 
shall use the averaged equations instead, which become 



H 2 — Q m0 e N 
P($' + $) + 



1 - e iN 
+ n r0 + -$ 2 + T^Ucff^) , C!2) 



3 4N dV e 



3H 2 



ft 



H 2 



. 



We now integrate the equations from N = — 70 to N = 
5, ranging from the Planck epoch to some time in the 
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FIG. 1. Evolution of log pdm for m — 10 -55 (dot-dashed), 
m — 10~ 35 (dashed) and m = 10 -15 (dotted), which is equiv- 
alent to the noninteracting compared to the back- 
ground density, log(p m + pr) (solid). 
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log(a) 

FIG. 3. Evolution of log pdm. for m — 10 -60 with the initial 
conditions \i — 1 (double-dot-dashed), \i — 2.609 (dashed) 
and x* — 10 (dot-dashed), as compared to background density 
log(p m + p r ) (solid) and to the noninteracting case (dotted). 
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log(a) 

FIG. 2. Evolution of log pa e for m = 10 -55 (dot-dashed), m = 
10 -35 (dashed) and m = 10 -15 (dotted), which is equivalent 
to the noninteracting scenario, as compared to background 
density, log(p m + p r ) (solid). 



future. We present our results in terms of log a instead 
of N, since it is easier to translate the former into redshift 
z. Nucleosynthesis takes place around log a = — 10. 

We begin by studying variations of the mass m and 
later of the exponent A. The parameters of the DE po- 
tential are kept fixed at A = 0.01, 4>q = 28.6 and A = 9.5. 
These are chosen so that there is a local minimum of 
the potential at the present epoch, while also complying 
with the bound for A from nucleosynthesis. As long as 
they verify these conditions, they can be changed with- 
out influencing the qualitative results. The parameters 
of the polynomial P{4>) are kept at_B = C = Z3 = lin 
all the plots presented, since no relevant change occurs 
when they are varied; in fact, even an increase by several 
orders of magnitude is similar to a change in to or A. 

We show in Figs. [T] and [2] the effect of the interaction 
on pdm and pde for the indicated values of the bare mass 
to. Notice that the result for to = 10~ 15 is equivalent 
to the noninteracting scenario and for this reason, it is 



shown for comparison in all figures. This is the range for 
which we can find initial conditions leading to suitable 
parameters for the Universe at present and after passing 
through a DM-dominated phase. We see that in these 
cases, pde is not significantly changed. Moreover, as men- 
tioned above, there is no effect close to the present for 
the studied range of masses. Nevertheless, the solution 
for the lowest mass (to = 10~ 55 ) is interesting in that 
scaling occurs for both DE and DM. In that case, the 
oscillations do not start until the matter-dominated era, 
which may be problematic for structure formation. 

For much smaller masses, such as the one predicted by 
Eq. (j25|) (to = 10~ 60 ), this is even more dramatic, as 
shown in Figs. |3] and HI Two classes of solutions exist 
with an accelerated expansion at present, both of them 
unrealistic: in the first case dark matter never dominates 
(dotted case, Fig. [3J, while in the second it dominates, 
but does not oscillate around the minimum of its poten- 
tial (dot-dashed case, Fig. |3J). In this second type of 
solution, oscillations never start and the x field instead 
slowly rolls down its potential, thus giving rise to an ac- 
celerated expansion. Furthermore, there is a special so- 
lution of the first type (dashed) that possesses transient 
accelerated behavior, although with the same problems 
already mentioned. 

The two regimes, to > 10 -55 and to = 10 -60 , have a 
frontier, at to = 5.9 x 10~ 57 , for which the solution still 
agrees with the observational constraints and presents a 
period of transient accelerated expansion, as shown in 
Figs. [5] and [6] Besides being transient, this solution 
presents other features: it clearly shows the onset of oscil- 
lations to be very close to the present- which, in principle, 
impairs structure formation; fldm is non-negligible during 
nucleosynthesis, as is VLd e , but their sum is small enough 
to comply with the bounds in Ref. [28|, VL r > 0.95. Also 
visible in Fig. [S]is the averaging of the oscillations around 
log a = —1.3. 

We now turn to the case A =^ A, which we study by 
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FIG. 4. Evolution of logp de for m = 10 with the initial 
conditions \i = 1 (dotted), Xi = 2.609 (dashed) and Xi — 10 
(dot-dashed), as compared to background density log(p m +p r ) 
(solid) and to the noninteracting case (superimposed with the 
Xi = X, since it is indistinguishable). 
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FIG. 5. Results for m = 5.9 x 10 showing the evolution of 
the relative densities fi<j e (solid), f2 dm (dashed), Sl r (dotted) 
and Q m (dot-dashed). 

fixing the mass at to — 10" 15 ~ 1 TeV. We obtain some- 
what similar results to those already found for masses 
to > 1CP 55 for a large range of A, as seen in Figs. [7] and 
[SI As before, there are no relevant effects at the present 
epoch, with the DM density changing only at early times. 
However, there is a special case for A = 2.8, where a 
solution with transient acceleration is found again, as 
shown in Figs. [S] and [TUJ This transient solution does 
not present the problem found previously in the other so- 
lution of this kind nor does it have an unnaturally small 
mass. The only slight anomaly is that fldm starts in- 
creasing around log(a) = —8, an effect of the interaction. 

The same transient result is also found for other pairs 
of (A, to) and their values are plotted in Fig. [TT] The 
expression found for the "transient line" is 

A = -0.1625 log(m) + 0.3706 . (33) 

This expression is rather similar to Eq. (|23l) . changing 
A to A, since the slope is —2 In 10/^(0) ~ — 21nlO/0o = 



FIG. 6. Results for m — 5.9 x 10 _J showing the evolution 
of the deceleration parameter q (dashed) and the DE EOS 
parameter Wde (solid). Also shown is the effect of the oscil- 
lations on the deceleration parameter before log a = —1.3; at 
that moment the oscillations are averaged and henceforth the 
evolution of the relevant quantities is obtained in terms of Eq. 
(32]). 
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log(a) 

FIG. 7. Evolution of logp dm for m = 10~ 15 for A = 9.5 
(non-interacting case, dotted), A = 6.5 (dashed), A = 4.5 
(dot-dashed) and A = 2.8 (double-dot-dashed), as compared 
to background density log(p m + p r ) (solid). 



—0.16102. This is expected, since in these transient so- 
lutions we anticipate the interaction to be relevant just 
until the present. As a consequence, it is not surpris- 
ing that the smaller the mass, the greater is the required 
value for A. 

However, not all of those (A, to) pairs are equally in- 
teresting, since in situations with mass smaller than 
10~ 15 , such as to = 10 -20 , the DM density parame- 
ter starts growing sooner and can exceed the mentioned 
limit from nucleosynthesis. This trend is not verified for 
the lowest of masses, as is attested by the solution for 
to = 5.9 x 10~ 57 . 

The explanation for this transient behavior lies in the 
interaction: the addition of pa-m to the effective potential 
raises the minimum of Vd e , thus allowing for the field to 
escape and continue to roll down the exponential. How- 
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FIG. 8. Evolution of logp d e for m = 1(T 15 for A = 9.5 (non- 
interacting case, dotted), A = 6.5 (dashed), A = 4.5 (dot- 
dashed) and A = 2.8 (double-dot-dashed), as compared to 
background density log(p m + p r ) (solid). 
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FIG. 9. Results for m = 1CT 15 and A = 2.8 showing the 
evolution of the relative densities Qde (solid), Qdm (dashed), 
Q r (dotted) and Q m (dot-dashed). 
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FIG. 10. Results for m = 1CT 15 and A = 2.8 showing the 
evolution of the deceleration parameter q (dashed) and the 
DE EOS parameter Wde (solid). 



FIG. 11. Plot of the parameter space singling out the line 
where transient acceleration solutions occur and the line 
where M 2 ((f>o) = Hq, the limit for oscillations. 



ever, for parameter values below the "transient line" of 
Fig. [Til w e obtain unrealistic scenarios; for those values 
the interaction is stronger, so either the solution presents 
a very low DM density at present or the accelerated phase 
is nonexistent, both of which are in conflict with obser- 
vational data. 

Albeit appealing, this transient scenario requires some 
degree of fine-tuning, since the precision required is ~ 
10~ 3 in A, using bounds on the values of fijo in compli- 
ance with [33| . The reason for that is that we need to 
"dissolve" the minimum for the field to run, but still al- 
low for cj) to slow down enough to create an accelerated 
expansion phase. Such solutions have been found previ- 
ously in the absence of interaction, for example, in the 
model of Ref. [29] , for a similar degree of fine-tuning in 
Vde- The present result shows that our model provides 
an alternative path to achieve such solutions, similarly 
to the two-field quintessence approach of Ref. [3J| . The 
advantage of the present proposal is that it requires tun- 
ing only of A or to, the initial condition of x fixed to 
reproduce the present dark matter density. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have studied a model with two cou- 
pled scalar fields, in which one plays the role of dark 
matter and the other of dark energy at late time. We 
find solutions that reproduce the current observational 
data and show that the interaction is irrelevant at the 
present epoch in the case A = A [cf. Eqs. (fT5)l and (|28|)] 
for masses larger than m = 5.9 x 10 -57 ~ 10 _29 eV. For 
this last case the interaction gives rise to a transient stage 
of acceleration and may be detectable, specially by its ef- 
fects on structure formation (see e.g. Refs. (g. l35l [3fi|). 
an issue that should be addressed in a future work. 

The study of the case A ^ A led us to find similar be- 
havior for most of the values of A, with an irrelevant con- 
tribution of the interaction at late time. The difference 
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appears only for specific values in a line of the parameter 
space of A and m. Those solutions lead to a transient 
stage of acceleration, not found in the other cases. This 
stage is compatible with observations and provides a way 
out to the accelerating regime-a useful property for defin- 
ing suitable asymptotic states free from future horizons 
in fundamental theories such as string theory (see Ref. 
[34j and references therein). The "transient line" is also 



the dividing line between solutions with and without an 
accelerated expansion. 
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